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Solving an inverse problem for a 
parabolic equation with a nonlocal 
boundary condition in the reproducing 
kernel space 


M. Mohammadi, R. Mokhtari* and F. T. Isfahani 


Abstract 


On the basis of a reproducing kernel space, an iterative algorithm for 
solving the inverse problem for heat equation with a nonlocal boundary con- 
dition is presented. The analytical solution in the reproducing kernel space 
is shown in a series form and the approximate solution un is constructed by 
truncating the series to n terms. The convergence of vn to the analytical 
solution is also proved. Results obtained by the proposed method imply that 
it can be considered as a simple and accurate method for solving such inverse 
problems. 


Keywords: Inverse problem; Parabolic equation; Nonlocal boundary condi- 
tions; Reproducing kernel space. 


1 Introduction 


The problem of finding the solution of partial differential equations with 
source control parameter has appeared increasingly in physical phenomena 
such as heat transfer, thermoelasticity, control theory, population dynamics, 
nuclear reactor dynamics, medical sciences, biochemistry, etc. [1, 2, 3]. The 
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parameter determination in a parabolic partial differential equation from the 
over-specified data plays a crucial role in applied mathematics and physics. 
This technique has been widely used to determine the unknown properties 
of a region by measuring a specified location in the domain. These unknown 
properties such as the conductivity medium are important to the physical 
process but usually can not be measured directly, or very expensive to be 
measured [2, 3]. In general, these problems are ill-posed. Therefore a variety 
of numerical techniques based on regularization, finite differences, finite ele- 
ment and finite volume methods are given to approximate solutions of such 
problems [2, 3, 4]. 

In recent years all kinds of boundary conditions and over-specified condi- 
tions arise in the inverse problems which make them more and more difficult 
to solve. The integral over-specified condition arises from many important 
applications in heat transfer, termoelasticity, control theory, life sciences, etc. 
Some different partial differential equations with nonlocal boundary and over 
specified conditions can be found in [5, 6, 7, 8, 9, 10]. 

The theory of reproducing kernels [11], was used for the first time at 

the beginning of the 20th century by S. Zaremba in his work on bound- 
ary value problems for harmonic and biharmonic functions. This theory 
has been successfully applied for solving a bunch of problems, see e.g. 
[12, 13, 14, 15, 16, 17, 18] and references cited therein. The book [19] provides 
excellent overviews of the existing reproducing kernel methods. 
In this paper, a new algorithm for determining unknown solution and un- 
known control parameter of the parabolic inverse problem with nonlocal 
boundary and integral over-specified conditions based on the reproducing 
kernel space, is presented. The advantages of the approach must lie in the 
following facts. The approximate solution converges uniformly to the analyt- 
ical solution. The method is mesh free, easily implemented and it needs no 
time discretization. Also we can evaluate the approximate solution v,,(2, t) 
for fixed n once, and use it over and over. 


The rest of the paper is organized as follows. In section 2 we describe the 
governing equation. Several reproducing kernel spaces are defined in Section 
3. The method implementation and convergence analysis are prepared in 
Section 4. Numerical results are presented in section 5. The last section is a 
brief conclusion. 


2 Governing equation 


Consider the inverse problem of determination a pair of functions {v,p} in 
the following parabolic equation 


Ov Ov 


ET jae = PU t Fat) (x,t)€Q=(0,1I)x(0,7] (1) 
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with the initial condition 
v(x, 0) = v(x), xz € [0,1] (2) 
nonlocal boundary conditions 
v(0,t) = v(1,t), vz(1,t) = 0, 0<t<T (3) 


and the integral over-specified condition 


re u(a,t)dx = E(t), t € [0,7] (4) 
0 


where f(z,t), (x), and E(t) are known functions. 


The existence, uniqueness, and continuous dependence of the solution 
upon the data for this problem are demonstrated in [20]. 


After taking integration from both sides of the equation (1) and using 
integral over-specified condition, we obtain 


E'(t) + vz (0,t) — fo f(a, tde 
E(t) : 


(5) 


p(t) = 


Then we have the following model problem 


, 1 
be _ oy _ PW tealO) fe MeDdey + #4), (2,8) € 2 = (0,1) x (0,71 


v(z,0) = v(x), v(0,t) =v(1,t), v2(1,t) =0. 


After homogenizing the initial condition, we have 


hi Qa 1 f(a,t)dx—E' (t 
a oa + fo an dy = EUG, t,t, We) (x,t) EQ= (0, 1) x eae 
u(z,0)=0, u(0,t)=u(l,t), ue(1,t) =9, 


where 


(u(a,t) + p(@)) (Ue (0, t) +! (0) 


F(x, t, 1, Ws) = 
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3 Reproducing kernel spaces 


Definition 1. Let H be a real Hilbert space of functions f :Q— R. Denote 
by (-,-) the inner product and let || - || = ./(-,-) be the induced norm in H. 
A real valued function K (x,y) :Q x Q-—> R is called a reproducing kernel of 
HA if the followings are satisfied: 


(i) K,(a) = K(z,y)€ H for ally €Q, 
(ii) f(y) = (f(x), Ky(z)) for all f € H and for ally €. 


Definition 2. A Hilbert space H of functions on a set Q is called a repro- 
ducing kernel Hilbert space if there exists a reproducing kernel K of H. 


Remark 1. The existence of the reproducing kernel of a Hilbert space H is 
due to the Riesz Representation Theorem. It is known that the reproducing 
kernel is unique. 


Now, we define some useful reproducing kernel spaces. The corresponding 
reproducing kernels can be found by the usual technique in many articles in 
literature (see [13]). 


Definition 3. Wo[0,1] = {u(x)|u(x), u(x), u(x) are absolutely continuous 
in [0,1], u(x) € L?[0,1],u(0) = u(1),u/(1) = 0}. The inner product and 
the norm in Wo[0, 1] are defined respectively by 


(U,U)w. = Sy u (0)v (0) + i u)(x)v) (x)dx, u,v € Wol0, 1], (7) 


lellw, _ (U, Uwe WE W,(0, 1]. 


The space Wo[0, 1] is a reproducing kernel space and its reproducing kernel 
function is called R, (a). 


Definition 4. W,[0,T] = {u(t)|u(t),u/(t) are absolutely continuous in 
(0, T],u’(t) € L?[0,T], u(0) = 0}. The inner product and the norm in 
W, [0,7] are defined respectively by 


T 
(i a) = dH (0)e ©) + if u(t)" (t)dt, u,v € W, [0,7], 


and 
el, = (Uw, we Wi[0, 7]. 


The space W,|0,7] is a reproducing kernel space and its reproducing 
kernel function r,(t) is given by 
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st+ 4?—18 t<-s, 
rs(t) -{ : 4 as 


st + St — 15° t>s. 


Definition 5. W2[0,1] = {u(x)|u(x),u’(x) are absolutely continuous in 
[0,1], w(x) € L?[0,1]}. The inner product and the norm in W2[0,1] are 
defined respectively by 


(u,U}w, = » u (0) (0) +f u'(a)u"(x)dx, u,v € W2[0, 1, 


lls = dwar 4 € W201. 


The space W2(0, 1] is a reproducing kernel space and its reproducing kernel 
function Q,(x) is given by 


Qy (x) = 


ee a<y, 
3 


2 
l+yr+ ¥a—- iy u>y. 


Definition 6. W3[0,T] = {u(t)|u(£) is absolutely continuous in [0, T], u’(t) € 
L?(0,T]}. The inner product and the norm in W3[0,T| are defined respectively 
by 


ie 
(U,V) w, = u(0)v(0) +f u'(t)u'(t)dt, u,v € W3[0, 7), 


and 
lellw = (u, Uwe? Wwe W3(0, T]. 


The space W3[0,T] is a reproducing kernel space and its reproducing 
kernel function q,(t) is given by 


1l+tt<s, 
wld ={it t>s. 


Definition 7. W(Q) = {u(z, t) eas is completely continuous in Q, jones € 


L?(Q), u(x, 0) = 0, u(0,t) = u(1,t), ue(1,t) = 0}. The inner product and the 


norm in W(Q) are defined respectively by 
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2 - 2 gi 2 gi 
Oo O Oo” 0 
(u,U)wy = >, Ee yi U0, 4) Ore yi 0, o| dt 


0 
T rye @ a 9 : 
+f | = Or u(z, t) an savtesd) dxdt, u,v € W(Q), 


I|24|| = (U,U) ws Ue W(Q). 


and 


Theorem 1. W(Q) is a reproducing kernel space and its reproducing kernel 
function is 


Ky,s)(@,0) = Ry(x)rs(), 
such that for any u(x,t) € W(Q), 


u(y, 8) = (u(x,t), Kiy,s)(2;t)) ws 


where Ry(x), rs(t) are the reproducing kernel functions of Wo[0,1] and 
W; [0,7], respectively. 


Proof. see [19]. 


Definition 8. wa Q) = {u(a, ty|5 2 is completely continuous in Q, pea € 


Ox 


L?2(Q)}. The inner product and the norm in W(Q) are defined respectively 
by 


ile 0); eh Fee Na pavl0.d) dt 
ie (OV ins 


ao 0 ae 
af r 7 [ew a u(a,t) a2 sute.d)| dxdt, u,v €W(Q), 


and 


lull. = (uu, wew). 


Ww) is a reproducing kernel space and its reproducing kernel function is 


Gy,s) (x,t) = Qy(@)qs (t). 
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4 The method implementation 


By defining the linear operator L : W(Q) > W(Q) as 


du OPu : Hee f(a, t)dx — E'(t) 


cee ae Pe EW a 


model problem (6) changes to the following problem 


Lu(a2,t) = F(2,t,u,uz), (x,t) €Q, 


(8) 
u(z,0) = 0, u(0,t) = u(1,t), Uz (1,t) = 0. 
Lemma 1. FL is a bounded linear operator. 
Proof. see [13]. 
Now, we choose a countable dense subset {(21,t1), (v2, t2),..., } in Q, and 


define 
oi(z, t) = G(w,,t,)(2; 8); wi(z, t) 3 L* $i (x,t), 


where L* is the adjoint operator of L. The orthonormal system {~;(z, t)}2, 


of W(Q) can be derived from Gram-Schmidt orthogonalization process of 


{Wi(x, t)}F2y as 
wi(z, t) = S- Bind (as t), 
ke1. 


where the orthogonal coefficients 3;, are given by 


1 os = 
eu? t=k=1, 
er oth, 
2, 
|rn-Sg 
Bik = i-1 ae 
Scag Byx 
as itk 


9 


i-1 : 
lwel2— So 2, 
j=l 


where 
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Jj 
= (= Ainteosn(e)) . 
(a,t)=(ai,ta) 


m=1 


Like in [13], we get the following theorems. 


Theorem 2. Suppose that {(xj,t;)}%, is dense in Q, then {w;(x,t)}%, is 


a complete system in W(Q) and w(x, t) = Lig ayy. (OO) gees 


Theorem 3. If {(xj,ti)}92, is dense in Q, then the analytical solution of (8) 
is 


u(x,t) = S- > Biz [F (tx, th, U(x, te), Oru(O, te))] Vila, t). (9) 


i=1 k=1 


By truncating the series in (9), we can obtain the approximate solution of 
(8). But, since the the series terms are not known, we need to construct an 
iterative method for obtaining the approximate solution. For this purpose, 
we choose nonnegative integer n and put the initial function uo(z,t) = 0. 
Then the approximate solution is defined by 


i=1 
where 
Bi = >> Bix F(x, te, tes (te, te), Otx1(0, te). (11) 


k=1 


On account of (26), the approximate solution p,,(t) can also be obtained by 


om E'(t) + Optn (0, t) 0 -f fle,tde be 


4.1 Convergence analysis 


The convergence of u, (x,t) can lead to that of p,(t), due to (26). So we 
only need to show that the approximate solution u,, (x,t) converges to the 
analytical solution u(z,t). At first, the following lemma is given. 
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Lemma 2. Assume that Un is a bounded sequence in W(Q), Un ala a, 


(fn, tn) > (y,8), asn > co. If F(x,t,u(a,t), ux(0,t)) ts continuous, then 
F(n, tn, Un-1(En, tn), OxUn—1(0, tn)) = Pgs §, uly, 8), 0,u(0, s)). 
Proof. Similar to proof of Lemma 2 in [13], we have 


|Un—1(an, tn) — u(y, s)| +0, as noo. 


Since 


ltn — 5] < Van — yl? + [tn — 5I?, 
if follows that 
(0, tr) —> (0, s). 
Thus in a same manner 
|OxUn—1(0, tn) — Ort(0,s)| 70, as noo. 
The continuation of F(ax,t, u(x), v(#)) implies that 


F (an, tn, Un—1(@n, tn), Oxtn—1(0,tn)) > F(y, s, u(y, s), Oxu(0,s)), as noo. 


Theorem 4. Suppose that un is a bounded sequence in W(Q) and (8) has a 
unique solution. If {(ai,ti)}%21 is dense in Q, then the n-term approximate 
solution u,(x,t) derived from the above method converges to the analytical 


solution u(x,t) of (8) in W(Q), such that 
u(z, t) = S- Bw (a, t), 
i=l 


where B; is given by (11). 


Proof. Similar to proof of Theorem 4 in [13], u,(a,t) converges to u(a,t) of 
the form 


t=1 
such that 
Lii(x1, ti) _ F (a1, ti, ui—1(21, ti), Ox —1(0, ti)). 


Since {(2x;,t;)}92, is dense in Q, for each (y, s) € Q, there exist a subsequence 
{@n,, tn; }721 such that 


(2nj1tn;) + (y,8) G00). 
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We know that Li(an,,tn;) = F(&n,,tn;,Un;—1(En,;; tn; ), Oxtin;—1(0, tn,))- 
Let j > co, by Lemma (2) and the continuity of F’', we have 


(Lu)(y, s) = Fly, 8, uly, s), 0,U(0, 8), 


which indicates that u(x,t) satisfies (8). 


Theorem 5. Under the conditions of Theorem 4, the approximate solution 
un(x,t) and its derivatives 0’! un(a,t), i = 0,1,2, j = 0,1, converge uni- 
formly to exact solution u(x,t) and its derivatives O't/u(2,t), i = 0,1,2, 
j = 0,1, respectively. 


Proof. 


\Oof ten (x, t) _ Ont u(x, t)| = ion” (Un(y; s) _ uly, s), K(c,t) (y, 5)) wl 
_ |(un(y, s) _ uly, 5), Ont? Kw) (y, 8) wl 
< Ot? Ket) (¥, 8)|lw llun(y, 8) — u(y, 8) lw 


BO cell Og |p, nm — CO. 


a 


5 Numerical experiments 


To test the accuracy of the proposed method, two examples are treated in 
this section. The results are compared with the exact solutions. 


Example 1. Consider problem (25)-(4) with 


aS) 


& 


(x) = 2+ cos(2r2), 
(t) =1+ e*, 
f(a,t) =1+4 4° e~* cos(2r2). 


It is easy to check that the exact solution is 


{v(x, t), p(t)} = fe-*(1 + cos(2rx)), —1}. 


Using our method, we choose 81 points in the region 9, and obtain the 
approximate solution vgi(x,t). We have listed approximate versus exact so- 
lutions, along with the relative errors at some nodal points at time T = + in 
Tables 1-2 and at time T = $ in Tables 3-4. Numerical results are in good 
agreement with the exact solutions. In Figs. 1-2, we display the exact and 
approximate solutions of v at times T = i and T = 4, respectively. In order 
to verify the convergence of the exact solution and its partial derivatives to 
the approximate solution and its partial derivatives, we depicted the relative 
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errors graphs of v, ve_ and Vee_ at time T = ; for different values of n in 
Figs. 3-5, respectively. The results show that the errors becomes smaller as 
n increases. 


Example 2. Consider problem (25)-(4) with 
y(x) =1+4 cos?(2r2), 
E(t) = se +1, 
f(a,t) = —8n7e' + 1627 cos? (2mx) — t — te’ cos?(2mx) — 1. 


The exact solution is 
{u(a, t), p(t)} = {1 + e cos?(2rx), 14 t}. 


Taking T = ; and choosing 81 and 144 points in the region Q, we have listed 
approximate versus exact solutions, along with the relative errors at some 
nodal points in Tables 5-6 and 7-8, respectively. Numerical results are in 
good agreement with the exact solutions and the accuracy of approximate 
solution is getting better as n increases. In Fig. 6, we display the exact and 
approximate solutions of v at time T = i Relative error distribution of uv 
at time T' = 4 is also given in Fig. 7a. It is clear that the numerical results 
are in good agreement with the exact solutions. Artificial errors 10~? were 
introduced into the right end and conditional condition. It can be seen from 


Fig. 7b that the error never affects the results of the method. 


Example 3. Consider problem (25)-(4) with 


p(x) = 1+ cos(272), 
E(t) = exp(—(27)*t), 
f(a, t) = (21)? cos(2r2) exp(—(2m)?t) + 2t(1 + cos(2max) exp(—(2m)?t + 10#?). 


The exact solution is given by 


{u(x,t), p(t)} = {(1 + cos(2ra) exp(—(2m)t), (27)? + 2t exp(10t7)}. 
Relative error distribution of v at time T = $ is given in Fig. 8a. It can be 
noted from Fig. 8a that our results are in better accuracy than the results 
in [20]. In order to demonstrate the stability of our algorithm, we shall give 
a perturbation « = 10~? to the right side function f(x,t) and over-specified 
condition E(t). The relative error distribution of v at time T = 5 depicted 
in Fig. 8b shows that the method is stable and gives excellent approximation 
to the solution. 
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6 Conclusion 


In this paper, the reproducing kernel Hilbert space method was applied suc- 
cessfully for solving an inverse problem for a parabolic equation with nonlocal 
boundary condition. Proposed method is shown to be of good convergence, 
simple in principle, easy to program and easy to treat the boundary condi- 
tions. It seems that the method can also be applied to higher dimensional 
inverse problems. We leave this to our further works. 
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Table 1: Relative errors of u(a,t) for Example 1; n = 81, T = + 


(x,t) exact Vapp Relative errors (2, t) Vesact Vapp Relative errors 
(1, og) ~—-2-998001 2.998100 3.315676E-05 Fit) 2.445035 2.443523. 6.182374E-04 
4.745) 1 1.000715 7.151300E-04 1,4) 2.692963 2.691261 — 6.321103E-04 
(1,745) 2.980099 2.980922 2.758099E-04 3,4) 1.818731 1.8169752 9.652627E-04 
(45-75) «2.636866 2.637919 + 3.995758E-03 7y%) 2.481098 2.479660 —_5.795225E-04 
(4, a5) 904837 1.909188 2.283755E-03 3,4) 1.409365 1.408556 —_-5.742847E-04 
(1, 45) 2.809675 2.810196 1.853823E-04 4,2) 1 1.002348  2.347800E-03 
(3,4) 1.002317 2.317300E-03 3,2) 1.152927 1.152638 —_2.505188E-04 
(3,4) 447420 1.448317 6.202361E-04 4,2) 1.400369 1.399123 8.895521E-04 
(3) 2.764994 2.764715 1.009424E-04 3.4) 1.148738 1.148205 4.639945E-04 
(3,4) 441248 1.441835 4.073198E-04 3,4) 1 1.002458 — 2.458400E-03 
(3,4) 846482 1.845673 4.378733E-04 1,4) 2.557601 2.552739 —-1.901260E-03 
Table 2: Relative errors of p(t) for Example 1; n = 81, T = + 
t Pexact Papp Relative errors t Dexact Papp Relative errors 
wm el -1.000064  6.453900E-05 x -1 -1.000684 — 6.838260E-04 
wel -1.000669  6.689750E-04 4 -l -0.999358 6.4141 28E-04 
wo «OC-l -1.000585 5.854630E-04 # +1 -0.9984953 1.504941 E-03 
ob -1 -1.000343  3.435430E-04 2-1 -0.997999 — 2.000439E-03 
z -1 -1.000717 _7.166060E-04 2-1 -0.997132 2.868264E-03 
Table 3: Relative errors of u(a,t) for Example 1; n = 81, T = 4 

(a, t) Vesact Vapp Relative errors (2, t) Vexact Vapp Relative errors 
(1, qi) ~—-2-998001 2.998135 4.480185E-05 £4) 2.445035 2.445058 —9.422361E-06 
(3.7) 1 1.000089 8.926900E-05 1,4) 2.692963 2.693930 3.587980E-04 
(1,qa5) 2980099 2.981268 3.922174E-04 3,4) 1.818731 1.816356 —1.305423E-03 
(sb, 75) 2.636866 2.640780 1.484401E-03 (45,4) 2.481098 2.485261 1.677789E-03 
(4,45) 904837 1.912180 3.854839E-03 2,4) 1.409365 1.407343 1.434693E-03 
(1,75) 2.809675 2.809855 6.428644E-05  ($,2) 1 1.002219 2.219280E-03 
(4,4) 1.004321 4.321410E-03 3,2) 1.152927 1.151739 —1.030525E-03 
(2,4) 447420 1.447886 3.219329E-04 4,) 1.400369 1.407778 5.291163E-03 
(1, 4) 2.764994 2.765560 2.047726E-04 3,4) 1.148738 1.146675 1.795169E-03 
(2,4) 441248 1.441568 2.220845E-04 4.4) 1 1.001966 1.966360E-03 
(3,4) 846482 1.845236 6.746695E-04 1,4) 2.557601 2.556291 5.124707E-04 
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Table 4: Relative errors of p(t) for Example 1; n = 81, T = 4 


+ 


Pexact Papp 


Relative errors 


t 


Pexact Papp 


Relative errors 


T00 -1 -0.999794  2.063958E-04 + -1 -0.977383 2.261733E-02 
to Co -0.997471  2.528598E-03 x -1 -0.979758 2.024151E-02 
ie -1 -0.997903 2.096714E-03 x -1 -0.975690 2.430969E-02 
a -1 -0.976715 2.328518E-02 3 -1 -0.971106 2.889415E-02 
$ -1 -0.976639 2.336095E-02 i -1 -0.964804 3.519553E-02 
Table 5: Relative errors of u(a,t) for Example 2; n = 81, T = + 
(x,t) Veeck Vapp Relative errors (2, t) Vexact Vapp Relative errors 
(1, zag) ©: 2001000 2.000890 —_5.523337E-05 7.4) 1.590680 1.591406 4.562790E-04 
(4,25) 2.012072 2.001377 5.315484E-03 (1,4) —-2.181360 2.164421 7.765664E-03 
(1,745) —- 2.010050 2.008924 5.601139E-04 32) 1 1.002636  2.636500E-03 
(qb.a5) 1.723344 1.715849 4.348790E-03 (45, 4) 1.799418 1.801355 —1.076466E-03 
(iq) 1 0.996456 3.544300E-03 (3,4) 1.305351 1.276493 _2.210723E-02 
(1, #) 2.105171 2.089205 7.583953E-03 4,3) 2.248849 2.171894 3.421967E-02 
(3,4) 2.117519 2.073065  2.099328E-02 3,2) 1.817382 1.751413 3.629875E-02 
(2,4) 1.279380 1.256824 1.762984E-02 (4,2) 1.312212 1.303871 6.356607E-03 
(1,4) —- 2.183148 2.116888 7.622889E-03 (2,1) 1.840405 1.768371 3.914031E-02 
(2,2) 1.283287 1.260171 1.801297E-02 (4,1) 2.284025 2.199602 3.696255E-02 
(3,4) 1 0.999525 4.751000E-04 1,4) 2.284025 2.266547 —7.652462E-03 
Table 6: Relative errors of p(t) for Example 2; n = 81, T = + 

t Dexact Papp Relative errors t  Pezxact Papp Relative errors 
Tia 1.001000 1.000778  2.216753E-04 3 1.125000 1.177465 4.663580E-02 
‘ian 1.012000 1.009609 2.363016E-03 + 1.166667 1.138380330  2.424571E-02 
<0 1.010000 1.007937 2.042468E-03 + 1.200000 1.157718084 3.523493E-02 
i 1.100000 1.141458 3.768888E-02 5 1.222222  1.170887564  4.200091E-02 
3 1.111111 1.156422 4.077993E-02 # 1.250000 1.211209275 —3.103258E-02 
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Table 7: Relative errors of u(a,t) for Example 2; n = 144, T = + 
(a, t) Veract Vapp Relative errors  (z, t) Vexnact Vapp Relative errors 
(1, qdhg) «2001000 2.000941 2.983507E-05 £4) 1.590680 1.595549 + 3.060781E-03 
(3. 7aa5) 2.012072 2.008807 1.622859E-03 (1,4) —-2.181360 2.180958 1.846293E-04 
(1,qa5) 2010050 2.009449 = 2.992512E-04 3.4) 1 1.013353 1.335274E-02 
(joi) 1.723344 1.722972 2.159116E-04 (7,5) 1.799418 1.810589 6.207763E-03 
(4 1 1.001805  1.805350E-03 2,4) 1.305351 1.309086. 2.861377E-03 
(1,qq) 2.105171 2.101083 1.941670E-03 (5, %) 2.248849 2.230536 8.143059E-03 
(4,4) 2.117519 2.102535 7.076129E-03 3,2) 1.817382 1.808653 4.803084E-03 
(2,4) 1.279380 1.276005 2.637737E-03 4,3) 1.312212 1.325054 —9.786674E-03 
(1, 4) 2.133148 2.129999 1.476139E-03 (3,4) 1.840405 1.833820 3.578477E-03 
(2,4) 1.283287 1.280545 2.136414E-03 4,4) 2.284025 2.266462 _7.689720E-03 
(3,4) 1 1.008052 8.052170E-03 1,4) 2.284025 2.293042 3.947540E-03 


Table 8: Relative errors of p(t) for Example 2; n = 144, T = 4 


4 
€. 


t Pevact Papp Relative errors t  Dexact Papp Relative errors 
70 1.001000 1.000929 7.116752E-05 + 1.125000 1.148678 2.104750E-02 
aa 1.012000 1.010731951 = 1.253013E-03 i 1.166667 1.145867427 1.782820E-02 
he 1.010000 = 1.008977553 =: 1.012324E-03 z 1.200000 =1.196111194 3.240672E-03 
40 1.100000 = 1.118037 1.639767E-02 2 1.222222 1.264996 3.499686E-02 
3 1.111111 1.131717 1.854514E-02 + 1.250000 =1.245162209 3.870233E-03 
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Figure 1: Exact and approximate solution of v for Example 1 at T = 4 
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Figure 2: Exact and approximate solution of v for Example 1 at T = 4 
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Figure 3: Relative errors graphs of v for Example 1 at time T = re a(n = 36), b(n = 64 


c(n = 81), d(n = 100) 
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Figure 4: Relative errors graphs of vz; for Example 1 at time T = + a(n = 36), 
b(n = 64), c(n = 81), d(n = 100) 
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Figure 5: Relative errors graphs of vert for Example 1 at time T = + a(n 36), 


b(n = 64), c(n = 81), d(n = 100) 
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Figure 6: Exact and approximate solution of v for Example 2 at T = + 
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Figure 7: Relative error graphs of v 


for Example 2 at time T = 33 a(without noise), 
b(with noise) 
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Figure 8: Relative error graphs of v for Example 3 at time T = 33 a(without noise), 
b(with noise) 


